Report for Model 6

— for both jointly and separately

2024-08-11

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.02
## LDevsum       1.02       1.05
## dh0           2.36       3.84
## dh1           5.65      11.83
## dh2           3.62       7.09
## dl0           1.62       2.49
## dl1           2.70       5.11
## dl2           2.89       5.27
## dl3           1.22       1.57
## 
## Multivariate psrf
## 
## 4.85

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1184.41236 22680.2827
DIC3 1139.45554 20662.2708
PWAIC 33.77812 948.6484
WAIC 1158.38855 21010.0069

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.02
## LDevsum       1.02       1.06
## dh0           4.37       7.72
## dh1           7.31      12.81
## dh2           3.59       6.71
## dl0           1.42       2.11
## dl1           1.42       2.07
## dl2           1.12       1.33
## dl3           2.24       3.69
## 
## Multivariate psrf
## 
## 7.27

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1377.23924 24118.439
DIC3 1243.55301 21367.683
PWAIC 89.58375 1524.755
WAIC 1311.70933 22092.162